  % Пример на рупорный раскрыв волновода с волной H10
clc
k = 2*pi/lmbd; % Волновое число

M=10;  % Половинное число шагов интегрирование по осям
Angls=360; delta=2*pi/Angls;    % Шагов по угловой координате и шаг для DN

teta=-pi:delta:pi; % Полярный угол 

aaxr=7; % ('Ширина раскрыва вдоль H (по Х:)');
axr=aaxr * lmbd;  

bbyr=3.5; % ('Ширина раскрыва вдоль E (по Y:)');
byr=bbyr * lmbd;  

phi=phi0./180 * pi; % Азимутальный угол 

  %Матрица, соответствующая диаграмме направленности в плоскости (Z0PHI)
Frup=ones(size(teta));
%Расчет матрицы множителя системы
[Frup]=d_raskrv(teta,phi,Angls, axr, byr, lmbd, M);
maxFrup=max(abs(Frup));  Frup=abs(Frup./maxFrup);

% Расчёт коэфицентов а2, а4
   n0_5=round(atan(tan(gama0)/2)/delta);
            f0_5=Frup(n0_5+Angls/2+1);

   n_1=round(gama0/delta);
            f_1=Frup(n_1+Angls/2+1);

a4=(12+4*f_1-16*f0_5)/3
a2=f_1-1-a4
